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Abstract - We present a study of the effects of magnetic fields on fluctuations 
and correlations in hadron resonance gas model. We find significant changes in 
the fluctuations of net baryon number, electric charge and strangeness. This 
is also reflected in various fluctuation ratios along the freezeout curve. 


Introduction. — Heavy-ion collisions (HIC) 
are investigated both theoretically and experi¬ 
mentally to understand the properties of nuclear 
matter at extreme conditions. One of the most 
important issues addressed in HIC is the possibil¬ 
ity for nuclear matter to undergo a phase transi¬ 
tions to quark matter. At low baryon density and 
high temperature nuclear matter is expected to 
smoothly cross over [I] to a quark gluon plasma 
(QGP) phase. Whereas, at high baryon density 
and low temperature the system is expected to 
have a first order phase transition HU. 

The study of the effect of magnetic field on the 
phase transition has become a subject of intense 
research for last few years. The phase transition 
in Quantum Chromodynamic (QCD) system is 
usually expected to occur around the QCD en¬ 
ergy scales A ~ 200 MeV. So one should be in¬ 
terested in the magnetic fields with strength B ~ 
(200 MeV)^ ^ 2 X 10^® G. Non-central relativistic 
HIC may create extremely strong magnetic field 
(~ ~ 10^®G) due to the relativistic motion 


of the charged particles. The magnetic field (R) 
may reach up to order of O.lm^, and 15m^ 
for SPS, RHIC and LHC energies respectively [5]. 

Magnetic fields can induce many interesting 
phenomena in QCD matter. For example the chi¬ 
ral magnetic effect, i.e., electric charge separation 
induced by chirality imbalance, along an exter¬ 
nal magnetic field, which also results in P and 
CP violation [5117]. On the other hand, magnetic 
catalysis [5] and inverse magnetic catalysis wm 
can affect the phase diagram of QCD matter. 

Lattice QCD studies show that in the presence 
of magnetic field the critical temperature may in¬ 
crease mm- Effects on hadron mass modifica¬ 
tion has also been reported [T3]. The presence 
of magnetic fields may increase the fluctuations 
and correlations as well as elliptic flow coef¬ 
ficient [15] of hot QCD matter. Anisotropic elec¬ 
tric conductivity of hadronic matter is an effect 
of strong magnetic field which has been found in 
lattice studies m- Such anisotropy in conduc¬ 
tivity should create an anisotropy in the dilepton 
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emission rate with respect to reaction plane and 
should be observed [TTlfTS] . The final observables 
in relativistic HIC may also depend on the evo¬ 
lution of the magnetic field. Depending on the 
electrical conductivity of the medium the mag¬ 
netic diffusion time r^ag may vary from 0.3 fm 
to 150 fm [TH[5D]. So for large diffusion time a 
constant magnetic field may be considered. 

In view of the importance of magnetic fields, 
we plan to study its effect on fluctuations and 
correlations of strongly interacting matter, us¬ 
ing the Hadron Resonance Gas (HRG) model [?, 
[niiiMiT]. Here the confined phase of QGD is 
modelled as a non-interacting gas of hadrons and 
resonances. Recently, the hadronic EOS for non¬ 
zero magnetic fields has been studied within the 
HRG model [IS]. On the other hand, study of 
fluctuations and correlations of conserved charges 
is a reliable way to study the phase transition. 
These quantities behave quite differently in the 
hadronic and quark phases P5MD] . Moreover, 
the fluctuations are expected to be enhanced near 
phase boundary and are related to the critical 
behaviour of strongly interacting matter EBISa- 
Furthermore, by studying event-by-event fluctu¬ 
ations [23] as a function of beam energy [3il[35] it 
would be possible to learn more about the QCD 
phase diagram. The paper is organised as follows. 
First we discuss the HRG model in the presence 
of a magnetic held. Thereafter we present our 
results and hnally we conclude. 

HRG model in presence of magnetic 

field. — The grand canonical partition function 
of HRG may be written as sum of all the partition 
functions Zi of each hadron i, where 

\nZi=±Vgi f -^-^ln[l (1) 

Here V is the volume of the system, gi is the 
degeneracy factor, T is the temperature, Ei = 
\/p^ + is the single-particle energy, rrii is the 
mass and pi = Bips + Sips + QiPQ is the chem¬ 
ical potential. In the last expression, Bi, Si, Qi, 
are respectively, the baryon number, strangeness 
and charge of the particle, p'S are corresponding 
chemical potentials. The upper and lower signs 
corresponds to fermions and bosons respectively. 


We have incorporated all the hadrons listed in 
the particle data book [22] up to mass of 3 GeV. 

For a constant magnetic field B^, the single 
particle energy levels for charged and neutral par¬ 
ticles are respectively given by [371139] 



where qi is the charge of the particle n is any 
positive integer corresponding to allowed Landau 
levels, Sz are the components of spin s in the di¬ 
rection of magnetic field, and ki is the anoma¬ 
lous magnetic moment. For a given s, there are 
2s -I- 1 possible values of s^. The gyromagnetic 
ratios are taken as gh = 2\qh/e\ for all charged 
hadrons [lasD]- Therefore the pressure for 
hadron may be written as 

for Qi = 0, (4) 

n'T 

n -^0 

ln[l ± for Qi 7 ^ 0 (5) 

where g' is the degeneracy other than spin. There 
are several simplifying assumptions that we are 
using here. Firstly we are considering a non¬ 
interacting HRG. Secondly we are considering 
point-like hadrons with their masses remaining 
unaffected by the ambient temperature, chem¬ 
ical potentials and magnetic fields. These are 
common assumptions that are regularly used to 
study the strongly interacting matter in the con¬ 
text of heavy-ion collisions as may be found in the 
Refs. [?,[1TJ[22H1I] and references therein. Simi¬ 
lar assumptions are also used by various lattice 


P-2 












Exploring effects of magnetic field on the Hadron Resonance Gas 




0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

T (GeV) 



1 0.08 0.1 0.12 0.14 0.16 0.18 

T (GeV) 


(b) 


(c) 



T (GeV) 


(d) 



0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

T (GeV) 


(e) 



0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

T (GeV) 


(f) 


Fig. 1: Variation of a-nd with T at fj, — 0 for baryon number, strangeness and electric charge. 


QCD groups to set the scales of their results for 
various thermodynamic quantities [mill]. This 
is true also in the studies of magnetised neutron 
stars as may be found in Refs. [5711551H5M5] . 

Also, apart from the above thermal contribu¬ 
tion, there is also a vacuum part which is in gen¬ 
eral divergent and needs regularisation and renor¬ 
malisation. However, theories with particles hav¬ 
ing spin > 1 are not renormalizable [46] . This 
is probably reflected through the fact that the 
renormalized contributions from particles with 
spin > 1 in the HRG vacuum is negative ESI- 
Since many of the thermodynamic quantities in¬ 
cluding fluctuations and correlations of conserved 
charges are unaffected by the vacuum part, we 
have neglected its contribution altogether. 

Results. — The order fluctuation and 
correlation are given respectively by the diagonal 
and off-diagonal components of the susceptibility 


X 


jk 

xy 




( 6 ) 


where n = j + k, Hx/y is chemical potential 
for conserved charge x/y, with x,y = B / S / Q. 

Figure [T] shows the variation of second order 
(x^) and fourth order (y^) susceptibilities of dif¬ 
ferent conserved charges with T at y = 0. At 
low r, dominating contribution to xb comes 
from protons and neutrons. With increase in T, 
other heavier baryons populate and contribute 
to Xb and hence susceptibilities increase with 
increase in temperature. Since all the baryons 
have baryon number ±1, there is no difference 
in magnitude for higher order susceptibility for 
Xb- Both Xb ^ad Xb increase with increase in B 
especially at high temperature. The number den¬ 
sity of protons increase with the magnetic field. 
At high temperature the A particles get excited. 
The number density of A+'''(1232) particles in¬ 
crease with magnetic field even faster than that 
for the proton for two reasons. Firstly, as the 
spin of A++(1232) is 3/2 the effective mass of 
A++(1232) decreases with magnetic field [55] . 
Secondly, for a spin 3/2 and doubly charged par¬ 
ticle the effect of degeneracy and magnetic field 
make the density increase faster. As a result the 
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Fig. 2; Xbs, Xqs, Xbq as a function of temperature at /r = 0. 


susceptibilities increase perceptibly at high tem¬ 
perature with the increase in magnetic field. 

The dominant contribution to xs at low tem¬ 
peratures come from kaons which have strange 
quantum number ±1. As a result magnitudes 
of Xs ^ad Xs similar in low T region up 
to 0.1 GeV. Since are spin zero particles, 
their populations get suppressed in presence of 
B. However, since T is small, population it¬ 
self is small and hence this suppression is also 
small. With increase in T, other strange mesons, 
like K* (spin one), starts populating the sys¬ 
tem. In presence of magnetic field populations 
of K*^ increases. As the contributions of 
and K*^ to the susceptibilities are opposite in 
nature, in presence of H, np to certain tempera¬ 
ture (T < 0.12GeV), susceptibilities for strange 
quantum number do not get affected much by the 
magnetic field. At higher temperatures, other 
charged strange hadrons (like start populat¬ 
ing the system. Furthermore, at high tempera¬ 
ture and high magnetic field, the contribution of 
iF*^, to the susceptibilities, is much more com¬ 
pared to that for K^. This also makes the sus¬ 
ceptibilities increase with magnetic field. Also, 
at high T, several strange baryons like S, are 
excited which have strange quantum number ±2 
and ±3 respectively. Therefore, in the high T re¬ 
gion, magnitude of xs increase for higher orders. 

In the low temperature domain, Xq and Xg 
decrease in presence of magnetic field. Beyond 
certain temperature, those increase with increase 
in magnetic field. Also this limiting temperature 


decreases with increase in magnetic field. The 
dominant contribution to the susceptibilities of 
conserved electric charges at low temperatures 
come from pions which are spin zero particles. 
Next dominating contributions are coming from 
kaons which are also spin zero particles. There¬ 
fore, the populations of these particles get sup¬ 
pressed in presence of magnetic held. As a re¬ 
sult susceptibilities of conserved electric charges 
decrease with increase of magnetic held at low 
temperature. With increase of temperature p, 
iF*, p, A etc., which are spin non-zero particles, 
appear in the system and these particles cause 
the susceptibilities to rise at high temperature. 

In hgure[2]we have plotted the correlations (oh 
diagonal susceptibilities) of conserved charges. 
For all the cases the correlations increase with 
magnetic held. The maximum ehect is observed 
in the baryon charge sector as the baryons with 
higher spin states and charge contribute in this 
sector. Here the leading contributors are p, 
and A. For all these particles the number density 
increases with increase in magnetic held, espe¬ 
cially at high temperature, as a result the corre¬ 
lation increases. The leading order contribution 
to BS comes from A and As the number 
densities of these particles increase in presence of 
a magnetic held the correlation also increases, es¬ 
pecially at high temperature. For the QS sector, 
the correlation is suppressed at low temperature, 
in presence of a magnetic held, as the major con¬ 
tributions come from . The number density of 
AT^, being a spin 0 particle, decreases in presence 
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Fig. 3: Products of moments of net-proton, net-charge and net-kaon as fnnctions of centre of mass energy. 


of a magnetic field. As the temperature increases 
the spin 1 particle K*^ appear in the medium. 
This makes the QS correlation increase at high 
temperature. 

Beam energy dependence of the products of 
moments. Experimentally measured moments 
such as mean (M), standard deviation (cr), skew¬ 
ness {S) and kurtosis (k) of conserved charges are 
used to characterize the shape of charge distri¬ 
bution. The products of moments can be linked 
with susceptibilities by the following relations, 


X 

X 


2 

X 
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X 
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X.X _ Q rr 
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— — tv 
2 ~ ^x^x' 

A.X 


(7) 


These ratios are independent of the volume of 
the system and play crucial role for the search 


of possible critical point in the QCD phase dia¬ 
gram. To make contact between our model and 
experimental data we need a parametrization of 
T, /r’s and B with ^/s. Such parainetrizaion ex¬ 
ist for multiplicities of identified hadrons in the 
absence of magnetic field [2^124] for central colli¬ 
sions. Here we use this parametrization to check 
effects of magnetic fields on the fluctuation ratios 
of net proton and net charge data given in |34I35) . 
However with this parametrization, a one-to-one 
comparison with experimental data is not ex¬ 
pected to be perfect as the system may not have 
equilibrated completely, leading to differences of 
fluctuations from those at thermal and chemical 
equilibrium. Secondly the parameters for central 
collisions are not same as peripheral collisions. 
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Lastly the effects of magnetic field has not been 
included in the above parametrization. We there¬ 
fore seek a qualitative comparison of the fluctua¬ 
tions with non-zero magnetic fields vis-a-vis those 
at -B = 0. A further example for possible future 
net-kaon fluctuation ratios is also presented. 

Figure [3] shows the of products of moments for 
net-proton, net-charge and net-kaon as functions 
of centre of mass energy ^/s. The two ratios 
{a'^fM)x and {Scr)^ involve odd derivatives xi 
and X 3 respectively. These should ideally vanish 
at zero densities. On the other hand, for large 
densities these approach 1. On the other hand, 
the ratio {Ka‘^)x is almost always close to 1 as 
obtained in HRG. This defines the general be¬ 
haviour of these ratios as functions of y/s. 

With non-zero magnetic field We find that 
/M)p is very sensitive with increasing ^/s. For 
(5 '(t)p the magnetic field effect is non-monotonic 
with y/s. Both these behaviours are prominently 
effects from the anomalous magnetic moment of 
the proton. We checked that without this anoma¬ 
lous part these ratios have insignificant depen¬ 
dence on the magnetic field. For {K,a^)p there is 
a small magnetic field dependence only at very 
low ^/s. For net-charge we find {a'^/M )q to have 
similar dependence on magnetic field as that of 
{a^/M)p, bnt the magnitude is now larger due 
to contributions of the low mass charged mesons. 
{Sa)Q is found to be most sensitive to the mag¬ 
netic field for lower -y/s. In these cases the effect 
of anomalous magnetic moment is subdominant. 
The dependence of (Kcr^)Q on magnetic field is 
small and of the same order as corresponding ra¬ 
tio of net charge. 

By construction, the 0-5% centrality data is 
close to the HRG results for B = 0. However the 
sensitivity of B underlines its importance to be 
considered as a freeze-out parameter in future. 
On one hand it seems that most magnetic field 
effects would be observed for peripheral heavy- 
ion collisions becanse large number of specta¬ 
tor particles would give rise to large magnetic 
field. On the other hand a larger system cre¬ 
ated in the central collisions may sustain even 
somewhat smaller magnetic field for a consider¬ 
able time m- In general it seems that the mag¬ 
netic field affects the lower order finctuation ra¬ 


tios the most. Therefore /M)x are suitable 
observables for fixing the magnetic field value in 
the HRG model. 

Finally we present a case for net-kaons. As we 
saw that the strangeness fluctuations are weakly 
dependent on the magnetic field, so is the case 
for all the ratios {a‘^/M)k, {Sa)k and {Ka^)k- So 
strangeness fluctuations are not much useful for 
detection of effects from magnetic field. We note 
here that no anomalous magnetic moments were 
introduced for the mesons as no conclusive ex¬ 
perimental data exist for them. 

Conclusion. — We have studied the fluc¬ 
tuations of conserved charges, namely baryon, 
strangeness and electric charge, using HRG 
model in presence of a magnetic field. We have 
chosen three values of magnetic fields, highest 
of which is close to the predicted magnetic field 
for non-central HIC at LHC energy. The study 
was done with several assumptions like neglect¬ 
ing interactions among the hadrons, ignoring the 
mass and size modifications in presence of tem¬ 
perature, chemical potential and magnetic field, 
as well as assuming that the effect of magnetic 
fields snrvive long enough for the hadrons to get 
affected by it. This study basically sets an upper 
limit for the magnetic field effects. 

The baryon number and electric charge num¬ 
ber susceptibilities are more sensitive to the mag¬ 
netic field as compared to the strangeness num¬ 
ber susceptibilities. Among the correlators, the 
baryon—charge correlation is found to be the 
most sensitive. The ratios of fluctuations that are 
measurable experimentally also manifest some ef¬ 
fects of the magnetic field. We loosely chose ther¬ 
modynamic parameters fitted for multiplicity at 
central collisions for B = 0 to see the deviation 
in an external magnetic field. We found that the 
fluctuation ratios for the net kaon have a soft de¬ 
pendence on the magnetic field whereas for net 
proton and net electric charge the dependence is 
quite strong. A reparametrization of freeze-out 
conditions with B as a parameter may be in or¬ 
der and will be stndied elsewhere. 

Acknowledgement. — We thank GSIR, 
DST and AvH foundation for support. SS thanks 
S. Das, S. Prasad, M. Younus and S. Maity for 


p-6 



Exploring effects of magnetic field on the Hadron Resonance Gas 


discussion. 


REFERENCES 

[1] Aoki Y., Endrodi G., Fodor Z., Katz S. D. 
and SzABO K. K., Nature, 443 (2006) 675. 

[2] Asakawa M. and Yazaki K., Nucl. Phys., 
A504 (1989) 668. 

[3] Ejiri S., Phys. Rev., D78 (2008) 074507. 

[4] Bowman E. S. and Kapusta J. L, Phys. Rev., 
C79 (2009) 015202. 

[5] Skokov V., Illarionov A. Y. and Toneev 
V., Int. Jr. Mod. Phys., A24 (2009) 5925. 

[6] Schafer T. and Shuryak E. V., Rev. Mod. 
Phys., 70 (1998) 323. 

[7] Fukushima K., Kharzeev D. E. and War- 
RINGA H. J., Phys. Rev., D78 (2008) 074033. 

[8] Shovkovy I. A., Lect. Notes Phys., 871 (2013) 
13. 

[9] Preis F., Rebhan A. and Schmitt A., JHEP, 
1103 (2011) 033. 

[10] Bruckmann F., Endrodi G. and Kovacs T. 
G., JHEP, 1304 (2013) 112. 

[11] D’Elia M., Mukherjee S. and Sanfilippo F., 
Phys. Rev., D82 (2010) 051501(R). 

[12] Bali G. S., Bruckmann F., Endrodi G., 
Fodor Z., Katz S. D., Krieg S., Sghafer 
A. and SzABO K. K., JHEP, 1202 (2012) 044. 

[13] E.V. Luschevskayaa, O.E. Solovjevaa, O.A. Ko¬ 
chetkova, O.V. Teryaev, Nucl. Phys. B896 
(2015) 627. 

[14] Fu W., Phys. Rev., D88 (2013) 014009. 

[15] Mohapatra R. K. , Saumia P. S. and Srivas- 
TAVA A. M., Mod.Phys.Lett, A26 (2011) 2477. 

[16] Buividovich P. V., Ghernodub M. N., 
Kharzeev D. E., Kalaydzhyan T., 
Luchevskaya E. V. and Polikarpov 
M. I., Phys. Rev. Lett, 105 (2010) 132001. 

[17] Bratkovskaya E. L., Teryaev O. V. and 
Toneev V. D., Phys. Lett, B348 (1995) 283. 

[18] Gupta S., Phys. Lett, B597 (2004) 57. 

[19] Tuchin K., Phys.Rev., C82 (2010) 034904. 

[20] Deng W. and Huang X., Phys. Rev., C85 
(2012) 044907. 

[21] Braun-Munzinger P., Stachel J., Wessels 
J. P. and Xu N. , Phys. Lett., B344 (1995) 43. 

[22] Cleymans j., Elliott D., Satz H. and 
Thews R. L., Z. Phys., C74 (1997) 319. 

[23] Cleymans J., Oeschler H., Redlich K. and 
Wheaton S., Phys.Rev., C73 (2006) 034905. 

[24] Karsgh F. and Redligh K., Phys. Lett, B695 
(2011) 136. 


[25] Endrodi G., .JHEP, 1304 (2013) 023. 

[26] Bhattacharyya A., Das S., Ghosh S. K., 
Ray R. and Samanta S., Phys. Rev., C 90 
(2014) 034909. 

[27] Bhattacharyya A., Ray R., Samanta S. and 
SUR S., Phys. Rev., C 91 (2015) 041901(R). 

[28] Jeon S. and Koch V., Phys. Rev. Lett., 85 
(2000) 2076. 

[29] Asakawa M., Heinz U. W. and Muller B., 
Phys. Rev. Lett., 85 (2000) 2072. 

[30] Ejiri S., Karsch F. and Redlich K., Phys. 
Lett, B 633 (2006) 275. 

[31] Stephanov M. a., Rajagopal K. and 

Shuryak E. V., Phys. Rev. Lett., 81 (1998) 
4816. 

[32] Hatta Y. and Ikeda T., Phys. Rev., D 67 
(2003) 014028. 

[33] Stephanov M. A., Rajagopal K. and 

Shuryak E. V., Phys. Rev., D 60 (1999) 
114028. 

[34] Adamczyk L. et al., Phys. Rev. Lett., 112 
(2014) 032302. 

[35] Adamczyk L. et al., Phys. Rev. Lett., 113 
(2014) 092301. 

[36] Olive K. A. et al.. Review of particle physics, 
C38 (2014) 090001. 

[37] Broderick A., Prakash M. and Lattimer 
J. M., Astrophys. J., 537 (2000) 351. 

[38] Broderick A. E., Prakash M. and Lat¬ 
timer J. M., Phys. Lett. B, 531 (2002) 167. 

[39] Bjorken j. D. and Drell S. D. , Relativistic 
Quantum Mechanics (McGraw-Hill: New York) 
1964. 

[40] Ferrara S., Porrati M. and Telegdi V. L., 
Phys.Rev., D46 (1992) 3529. 

[41] Bazavov a. et Ah., Phys. Rev., D 86 (2012) 
034509. 

[42] Borsanyi S. et Ah., JHEP, 1011 (2010) 077. 

[43] Mallick R., Schramm S.,Dexheimer V. and 
Bhattacharyya A., Mon. Not. Roy. Astron. 
Soc., 449 (2015) no.2 1347. 

[44] Ghakrabarty S., Bandyopadhyay D., and 
Pal S., Phys. Rev. Lett, 78 (1997) 2898. 

[45] Bandyopadhyay D., Ghakrabarty S. and 
Pal S., Phys. Rev. Lett, 79 (1997) 2176. 

[46] Peskin M. E. and Schroeder D. V., An In¬ 
troduction to quantum field theory (Addison- 
Wesley, Reading) 1995, p. 842. 


P-7 



